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Various factors can be crucial to hydrodynamic instabihty.: shear, gravity, surface tension, 
heat, centrifugal force, etc. Some instabilities lead to a different flow; some to turbulence. 

We only discuss the linearized analysis of instability of a few sample problems. More 
examples will be discussed in later chapters. 

Let us take the parallel flow as the example and outline the standard procedure of lin- 
earized analysis as follows: 

1. Assume the flow domain to be of inflnite extent in the x direction. Let the basic flow 
have velocity U {z) in the x direction. 

2. Derive the linearized equation for infinitestimal perturbations and homogeneous bound- 
ary conditions. 

3. Consider sinusoidal wave- like disturbances so that all unknowns are of the form 

fi^fiiy,z) e'^'^-^'^ 

4. Deduce the partial differential equation(s) in the transverse plane {y, z) and obtain an 
eignevalue problem. 

5. Solve for the eigenvalue cu for given real k. 

6. If the eigenvalue is complex, cu — cUr + iuJi, find the condition under which a;^ > 0. Since 
the flow is unstable a;^ > as time grows, stable if cui < 0, and neutrally stable if 

LJi = 0. 
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5.1 Kelvin-Helmholtz instability of flow with discon- 
tinuous shear and stratification. 

5.1.1 Nonlinear formulation 

Consider two immiscible fluids. The upper {z > 0) and lower (z < 0) fluids have different 
densities pi and p2- Before any perturbations appear the lower fluid is stationary, while the 
upper fluid moves at the steady and uniform velocity (jo = Ui. With gravity, this is called 
the Kelvin-Helmholtz instability problem. For the time being, let us assume that there is no 
gravity. 

Now let there be a disturbance on the interface whose vertical displacement is 2; = C(^) 
The corresponding disturbances of velocity and pressure are qi and pi in the upper fluid and 
g*2 and p2 in the lower fluid. 

Let us flrst formulate the nonlinear problem. In the upper fluid, we have 

^i = V$i, and V^$i=0, z>0 (5.1.1) 
and the Bernoulli equation, 

^|l,(V^ = _Pl + C. (5.1.2) 
at 2 pi 

In the lower fluid, we have : 

= V$2, and V^$2 = z<0 (5.1.3) 

with the Bernoulli equation, 

?ll^^^ = -'^ + C, (5.1.4) 

at 2 P2 

On the interface F — z — (^{x, t) — the kinematic conditions are 

dF 

— + V$i-VF = (5.1.5) 
at 



and 



or, more explicitly. 



dF 

— + V$2-VF = (5.1.6) 
at 



dt dx dx dz 
dt dx dx dz ^ 



The dynamic condition is 

Pi = P2, or z = (. (5.1.9) 
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+ -C,J=p.(^^ + ^— or. = C. (5.1.10) 
Far away from the interface, 

V$i ^ C/i?, z ^ oo, (5.1.11) 

V$2 ^ ^ -oo. (5.1.12) 

_ — * — * 

Note that in the absence of any transient disturbances, V$i = Uii, V$2 = U2i, so that 



5.1.2 Linearization 

We now hmit ourselves to infinitesimal disturbances. SpecificUy, the interface displacement 
is infinitesimal when compared to other length scales. 

$i = C/ia; + 0i, and ^2 ^ U2X + (f)2, (5.1.14) 

with 0„ <S n = 1, 2 and linearize all the conditions. We then get 

V^$i=0, z>0, (5.1.15) 

nad 

V2$2 = 0, z < 0, (5.1.16) 

Note that 

fix, C, t) = fix, 0, t) + C^U=o + ■■■, (5.1.17) 
On the inteface we have ther kinematic conditions 

and the dynamic condition 

(d(j)i d(j)i\ (d^2.jjd^2\ ^ ^c: 1 

where (5.1.13) has been used. 
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5.1.3 Normal mode analysis 

Assume a sinusoidal disturbance 



(5.1.21) 



We follow the usual convention that the real part of the final solution is to be taken to 
represent physical quaitities. To satisfy Laplace equation we take 



^1 Q-^^ Q^i^^-^^) 



b2 e'^^e*^^^"'^*^ 



(5.1.22) 
(5.1.23) 



Here A, 0i, 02 ^-re unknown constants. 

Substituting into boundary conditions on the interface eqs (5.1.18-5.1.20) we get three 
homogeneous alegebric equations for A, 0i, 02, 



{—iuj + ikUi)A = —k4>i 

{—iuj + ikU2)A — A;02 
Pi{—ioj + ikUi)(l)i = p2{—iuj + ikU2)(t>2 



(5.1.24) 
(5.1.25) 



which can be written in matrix form 



-iu) + ikUi 
-iuj + ikU2 




k 









-k 




iuj + ikU2) 





In order to get non trivial solutions the coefficient determinant of the algebraic equations 
must vanish, leading to the eigenvalue condition 

pi{-cu + kUif + P2{-UJ + kU2f = 



or 



The solution is 



(pi + P2)uJ^ - 2k{piUi + p2U2)co + k'^ipiUl + Up2Ul) = 



_ j^ PlUl + P2U2 , .,^/Plp2\Ul - U2\ 



pi + P2 



pi + P2 



(5.1.26) 



The lower sign corresponds to instability. Thus a velocity discontinuity is always unstable. 
This is the limiting case of the continuous velocity profile with a point of inflection. 
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5.1.4 Physical interpretation 

Let us use the mathematical results to see the physics ( Batchelor, p 516fF). For simlicity 
we choose pi = P2 = P and Ui = —U2 = C/ > so that (5.1.26) reduces to 



±ikU 



(5.1.27) 



and is purely imaginary. The positive root (upper sign) signifies instability. 

The basic flow has a vortex sheet at 2; = 0. Due to the disturbance, the vorticity 
distribution is changed. Across the interface the total vorticity is approximately 



du' 



0+)-«'2(^ = 0- 



(901 (902 I 



dx 



dx 



Note that a positive vorticity vector points into the paper (as does the y axis); the circulation 
is clockwise in the x — z plane. Prom the solution (5.1.24) and (5.1.25) we get 



901 
dx 



ikx—iuit 



—i{—iuj + ikUijAe 



ikx—iuit 



—i{—iuj + ikU\)C, 



and 



dx 



lo = ik(t)2e'''''-'^^ = i{-iuj + ikU2)Ae'''''-'''^ - i{-iuj + ikU2)C, 
Therefore the vorticity along the interface is 

[-2uj + k{Ui + U2)\C, = T2ikUC = 2kUAexp{ikx ^ i7i/2) ± kUt) 



Consider only the upper signs for the unstable solution, 



-2ikUC = 2kUAexp 



ik i X 



71 

2k 



= 2A;C/Aexp 



ik X 



+ kUt 



+ kUt 

since k — 2'k/X where A is th wave length. In real form the vorticity disturbance varies as 

2kUA cos 



k[x-\ 



^kUt 



whereas the interface varies as 



A cos kx e^^* 



The vorticity disturbance varies along the interface and leads the interface displacement C, 
in phase by 7r/2 or one quarter of a wavelength. As sketched in the upper part of Figure 
5.1.1, the vorticity disturbance is at the positive (clockwise) maximum at points like A where 
= and < and at the negative (counterclockwise) maximum at points like B where 
= and C^x > ^1 being zero at the crests and troughs of the interface. All these vortices 
tend to lift the interface crests and supress the troughs. 
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Figure 5.1.1: Disturbance vorticity along a velocity dicontinuity. Ui — U > 0,1/2 — —U < 0. 
Moreover the disturbances are convected by the average velocity 




= k{Ui - U2)C = 2UC 



which is positive near the crests and negative near the troughs. Therefore clockwise vorticity 
disturbances accumulate around A, while the counterclockwise disturbances are swept away 
from and thinned out around B. Because of this assymetry the wavy interface not only tend 
to amplify but to curl up in the forward direction as sketched in the lower part of Figure 
5.1.1. 

After some time the disturbance grows so large that nonlinearity must be accounted for. 
Numerical nonlinear analysis as well as laboratory experiments show that the interface curls 
up. 
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Homework K-H instability: Show that when gravity and surface tension are both 
present, the eigenvalue condition is 

_ PiUi + P2U2 

Discuss the condition for instability. 

Remark: The special case of f/i = f/2 = but g ^ 0, i.e., the instability of normally 
accelerated density discontinuity, is called the Rayleigh- Taylor instability. It is of interest 
when an interface is acelerated in the normal direction, and has a wide range of industrial 
and astrophysical applcations. 

Remcirk: The K-H instabilty problem can be generalized in a number of ways, by 
including viscosity and continuous density and/or velocity profiles, see Chandrasekhar. 



Plp2iUi - U2Y , 9 P2- Pi 



Tk 



(pi + P2y 



kpi+ P2 Pl+ P2 



(5.1.28) 



